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ABSTRACT 

We observed five central stars of planetary nebulae (CSPN) in the Large Magellanic Cloud (LMC) and 
three in the Small Magellanic Cloud (SMC) with the Far Ultraviolet Spectroscopic Explorer (FUSE), in the 
range 905 — 1 187 A. We performed a model -based analysis of these spectra in conjunction with Hubble Space 
Telescope (HST) spectra in the UV and optical range to determine stellar and nebular parameters. The signature 
of hot (T > 2000 K) circumstellar molecular hydrogen is found in the FUSE spectra of most objects. We also 
find evidence of X-rays in the wind of LMC-SMP 76. 

Subject headings: stars: AGB and Post-AGB — stars: atmospheres — stars: individual (LMC-SMP 1, LMC- 
SMP 29, LMC-SMP 50, LMC-SMP 76, LMC-SMP 83, SMC-SMP 1, SMC-SMP 3, SMC- 
SMP 5, SMC-SMP 22) — ultraviolet: stars 



1. INTRODUCTION 

Central Stars of Planetary Nebulae (CSPN) residing in the 
Magellanic Clouds (MCs) are an important asset to the study 
of CSPN evolution for two main reasons: their distances are 
well constrained, allowing their physical parameters (e.g., ra- 
dius and luminosity) to be determined with much less uncer- 
tainty than their Galactic counterparts, and the metallicities of 
their host galaxies are lower than in the Milky Way. Metal- 
licity is expected to influence a star's evolution by setting 
the efficiency of radiative driving during the phases when the 
star has a stellar wind: the asymptotic giant branch (AGB) 
phase, the post-AGB phase, and the [Wolf-Rayet] phase 
([WR]). Higher metallicity, and thus wind radiative acceler- 
ation, should increase the star's mass-loss rate and wind ve- 
locity, altering the layers exposed and chemical yields as well 
as the kinematics of the nebular shell. This has implications 
for galaxy evolution, through chemical enrichment and the 
dynamic interactions between the star's ejected material and 
the surrounding interstellar medium (ISM). The large fraction 
of stars that pass through these phases make their contribu- 
tion to the galactic chemical evolution significant (see, e.g., 
Marigo e t al.ll200 T] for a discussion). 

The far-UV wavelength range is useful for studying hot PN 
nuclei, as they emit the bulk of their observable flux in this 
range, and it contains many strong lines which are diagnos- 
tics of stellar paramet ers. It is also uncontaminated by nebular 
contin uum emission (Bi anchi et alJ Il997t [Herald & Bianchil 
l2004al) , which affects wavelengths longer than Lya, allowing 
us to characterize the physical parameters of both the nebu- 
lar and central star components through detailed quantitative 
modeling. 

Motivated by the above considerations, we 
dHerald & Bianchil l2004al — hereafter, HB04) performed 
a spectroscopic analysis on far-UV and UV data of seven 
CSPN in the Large Magellanic Cloud (LMC) using data taken 
with the Far-Ultraviolet Spectroscopic Explorer (FUSE). 
As a continuation of that work, we have obtained FUSE 
observations of five additional LMC CSPN (Bianchi's D034 
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DATA. FUSE IS OPERATED FOR NASA BY THE JOHNS HOPKINS UNI- 
VERSITY UNDER NASA CONTRACT NAS5-32985. 



program) and three Small Magellanic Cloud (SMC) CSPN 
(Bianchi's C056 program). These targets were selected 
from those expected to be the brig htest (in the far-UV) from 
the ph otoionization studies of [Domta & Meatherinaham 
dl991allbh : iDopita et alJ d!997l) : IVassiliadis et alJ dl996l 
1998a). The additional LMC targets (D034) were selected to 
expand the temperature coverage of the HB04 LMC sample 
to higher temperatures, and with the SMC sample we gain 
objects in a much lower metallicity environment. Together 
with the LMC sample of Bianchi's FUSE proposal B001, 
analyzed by HB04, these programs cover all the known 
CSPN in the MC observable by FUSE. Our analysis is similar 
to that of HB04, and the reader is referred to that work for 
more details and background information. 

This paper is arranged as follows. The observations and 
data reduction are described in § [2] A comparison of the 
spectra of the objects is presented in § [3] Our models and pa- 
rameter determinations are described in §|U Individual object 
results are presented in § [5] The implications of our results 
are discussed in §[6] and our conclusions in §|7] 

2. OBSERVATIONS AND REDUCTION 

For this work, we used far-UV spectra taken with FUSE and 
archive UV and optical spectra taken with the Hubble Space 
Telescope 's (HST) Faint Object Spectrograph (FOS) or Space 
Telescope Imaging Spectrograph (STIS), which we now de- 
scribe. 

2.1. Far-UV Spectra 

The FUSE observations (programs C056 and D034) of our 
sample stars are summarized in Table Q] They are some of 
the dimmest stellar objects yet observed by FUSE, and neces- 
sitated long integration times. FUSE covers the wavelength 
range 905 — 1 187 A at a spectral resolutio n of R ss 20 000 ( 
15 kms" 1 ). The instrument is described by Moos et alJ d20Qi 
and it s on-orbit performance is discussed by Sahnowet 
(2000). FUSE collects light concurrently in four different 
channels (LiFl, LiF2, SiCl, and SiC2), each of which is di- 
vided into two segments (A & B) recorded by two detectors, 
covering different subsets of the above range with some over- 
lap. We used FUSE'S LWRS (30"x30") aperture. These 
data, taken in "time-tag" mode, have been calibrated using 
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the FUSE data reduction pipeline, efficiency curves and wave- 
length solutions (CALFUSE v2.4). 

The spectra of these targets suffer from significant contami- 
nation from terrestrial airglow lines (such as O I and N I), and 
scattered solar emission features (O VI AA1032,38 and C III 
Al 175), which often interfere with the spectroscopic analysis 
by obscuring the stellar features. For the LMC sample, we use 
"night-only" data (i.e., taken when the spacecraft was on the 
dark side of the Earth), which minimize this contamination, 
although some residual airglow lines (typically O I) can still 
appear (an example of the difference between day and night 
observations can be found in HB04). 

Because of their faintness, most of our targets required 
many orbit-long exposures, each of which typically had low 
count-rates and thus signal-to-noise ratios. Each calibrated 
extracted sequence was checked for unacceptable count-rate 
variations (a sign of detector drift), incorrect extraction win- 
dows, and other anomalies. Segments with problems were not 
included in further steps. The default FUSE pipeline attempts 
to model the background measured off-target and subtracts it 
from the target spectrum. In some cases, the background is 
taken to be the model scattered-light scaled by the exposure 
time. 

All the "good" calibrated segments/exposures were com- 
bined using FUSE pipeline routine and overlapping ranges of 
the different channels were then compared for consistency. 
There are several regions of overlap between the different 
channels, and these were used to ensure that the continuum 
matched. For the fainter targets, often the continuum levels 
for different detectors did not agree. The procedure we typi- 
cally adopted was to trust the LiFla absolute flux calibration 
(which is the most reliable — Sah now et al.ll2.QO0t) , and scale 
the flux of the other detectors if needed. The LiFlb segment 
is sometimes affected by an artifact known as "the worm" 
(FUSE Data Handbook v 1 . 3) and was omitted if the effect was 
found to be severe. The SiCl detector seemed frequently off- 
target, and other segments offered higher S/N data in the same 
range, so its data was often not used. In the end, we typically 
used the SiC2 data for wavelengths below ~ 1000 A. Long- 
wards of 1080 A LiF2a data was employed. For the inter- 
mediate region (i.e., 1000-1080 A), data from LiFlb, LiF2b, 
or both were used (data from one segment was omitted if it 
was discrepant with that of the SiC detectors). The region be- 
tween 1083 and 1087 A is not covered by the LiF detectors, 
and as the SiC detectors in this range were off-target, we have 
omitted this region (appearing as gaps in the plots). Table Q] 
shows which segments and portions were used for the com- 
bined spectrum. Data from near the detector edges were also 
omitted if they looked inconsistent. 

Some individual objects warrant more discussion. The 
SMC observations had numerous difficulties. These obser- 
vations occurred during a period of time after FUSE lost one 
of its reaction wheels, affecting its pointing, but before jit- 
ter correction had been introduced into the FUSE calibration 
pipeline to account for such effects. For SMC-SMP 1, it ap- 
pears only the LiFl was aligned with the central star with any 
reliability (the fine error sensor camera, used for guiding, is 
part of the LiFl channel, so the location of the source is al- 
ways aligned in the LiFl apertures). Other detectors seemed 
to have gathered nebular emission lines with little or no stel- 
lar continuum for this object. We therefore only present the 
LiFl data for this object but note that the LiFlb (longer wave- 
length) segment is often affected by "the worm" (see above). 



Although the worm is not obviously present in the LiFlb ob- 
servations of this object, the low count rate prevents a defini- 
tive statement regarding this. During the observation of SMC- 
SMP 22, there was significant target drift as evidenced by in- 
constant count rates, and the LiF2 and SiC2 detectors seem 
to have been off-target. It also appears that another star was 
present in the aperture during much of the observation (dis- 
cussed in § 15.21 ). We present here data from the SiCl and 
LiFl detectors, but we stress that the absolute fluxes are un- 
certain due to what appears to be the target drifting in and out 
of the aperture throughout the observation. For SMC-SMP 5, 
no useful data was acquired. In the case of LMC-SMP 29, 
voltage and alignment problems resulted in no reliable FUSE 
data being collected. 

The FUSE spectra presented in this paper were obtained by 
combining the good data from different segments, weighted 
by detector sensitivity, and rebinning to a uniform dispersion 
of 0.05 A. 

2.2. UV and Optical Data 

Archival data from the HST's Faint Object Spectrograph 
(FOS) and Space Telescope Imaging Spectrograph (STIS) in- 
struments were also used longwards of 1200 A, as summa- 
rized in Table [2] The FOS archive data were taken in the 
1" aperture, whose actual diameter following the COSTAR 
installation is 0.86". The angular sizes of the LMC PN are 
typically <1", so the FOS spectroscopy includes the entire 
nebula (the exception is LMC-SMP 83, see Table [3]). The 
utilized FOS dispersers include G130H ( AA ~ 1 A), G190H 

( AA ~ 1 .5 A), G270H ( AA - 2~k) and G400H ( AA - 3 A). 
For LMC-SMP 83, medium-resolution STIS data is avail- 
able, taken through the 52"x0.2" aperture with the G140L 
( AA ~ 1 .3 A), G230L ( AA - 1 .5 A), G430L ( AA - 1 .6 A) 
and G750L AA - 1 .6 A) gratings. 

Generally, the flux levels of the FUSE and FOS/STIS data 
in the region of overlap are in good agreement for the LMC 
observations, except for LMC-SMP 83, which is known to be 
variable (§0). The SMC sample contains more discrepancies, 
however. For SMC-SMP 22, the UV flux level is much higher 
than that of the far-UV (discussed more in § 15.21) . For SMC- 
SMP 3 and SMC-SMP 1, the FUSE flux levels are higher by 
about a factor of 2 and 1.5 (respectively) than the FOS data. 
International Ultraviolet Explorer (IUE) spectra of these ob- 
jects display similar flux levels as those of FOS, therefore we 
we have scaled the FUSE flux levels to that of the FOS in 
these instances. 

3. DESCRIPTION OF THE SPECTRA 

The reduced FUSE spectra for the LMC sample, along with 
line identifications, are shown in Fig.Q] and the corresponding 
UV spectra are shown in Fig. 

Wind features are definitely seen in the far-UV spectra of 
most stars: S VI AA933,44, C III A977, and S VI AA933,44 
are visible in LMC-SMP 1 and LMC-SMP 76. The latter also 
shows C III A1175 and P V AA1 118,28. For LMC-SMP 50, 
the only feature obviously attributable to a stellar wind is 
O VI AA1032,38. LMC-SMP 83 shows also nebular emis- 
sion in O VI, implying a very hot central star (temperatures 
of T e ff > 125 kK are needed to gen erate these lines through 
photoionization — IChu et ai1l2004l) . In the far-UV spectra, 
H 2 abs orption makes the continuum level very difficult to set 

(§ED- 

In the UV range, LMC-SMP 50 shows some nebular fea- 
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tures and an underlying continuum, but strong stellar fea- 
tures are absent. LMC-SMP 83, LMC-SMP 1, and LMC- 
SMP 76 all display a P-Cygni profile of N V AA1238,43, 
as well as C IV AA1548,51 of different strengths. The blue 
edges of the P-Cygni profiles in both the far-UV and UV 
spectra indicate the terminal wind velocities of all objects are 
Voo ~ 500- 1000 kms" 1 , except in the case of LMC-SMP 83 
(Voo ~ 1600 kms -1 ) and SMC-SMP 3 ( > 2000 kms" 1 ). 

The reduced FUSE, and archive UV, spectra for the SMC 
sample are shown in Fig. [3] and Fig. |4] respectively. In the 
FUSE observations, nebular and airglow features contaminate 
the spectra (especially at short wavelengths). C III Al 175 is 
present in absorption in the spectra of SMC-SMP 1 and SMC- 
SMP 22. O VI AA1032,38 nebular lines are also seen in the 
latter. The only wind features in the UV spectra are C IV 
AA1548,51 seen in SMC-SMP 1. SMC-SMP 22 shows a more 
highly-ionized nebular spectrum than that of LMC-SMP 83. 

4. MODELING 

Modeling the central stars of the MC CSPN presents some 
challenges. First, their nebulae are typically compact (< 1"), 
and in most cases entirely contained in the FOS aperture used 
in the archival data available. Their UV and optical spectra are 
frequently contaminated by the nebular continuum and lines, 
which obscure and sometimes entirely mask the stellar spec- 
trum in these regions. Thus, one must rely on a smaller set of 
stellar features in the far-UV and UV to determine the stellar 
parameters. Second, the far-UV region, while not influenced 
by nebular continuum emission, is affected by absorption by 
electronic transitions of sight-line molecular hydrogen (Hi). 
In our objects, such absorptions cover the entire far-UV re- 
gion. The analysis of the spectra consists of modeling the 
central star spectrum, the nebular continuum emission long- 
wards of 1200 A, and the sight-line Hydrogen (atomic and 
molecular) shortwards concurrently and to find a consistent 
solution. However, we discuss each in turn for clarity. The 
reddening has to be derived concurrently, but it is small in all 
cases, and therefore it does not contribute much to the uncer- 
tainties. 

Through out this pap er, we adopt a LMC distance of D = 
50.6 koc dFeastl 119911) . an SMC distance of D = 60.0 kpc 
([Harries et all 120031) . an LMC metallicity of z = 0.4 Z Q 
( bufourll 19841) , and an SMC met allicity of z = 0.1 Z^, with 
values of "solar" abundances from lGrevesse & Sauvall(ll998l) . 

As wit h HB04, we a ssume Galactic foreground extinctions 
of 0.05 dBessell|[l991h . with t he remainder bein g from the 
Magellanic Cloud. We use the Fitzpatrick (1999) reddening 
law for the LMC and the lCardelli et al.l (1 19891) reddening law 
for the Galaxy (R v = 3.1) extinction, extrapolated down to far- 
UV wavelengths. 

4. 1. Molecular and Atomic Hydrogen 

Absorption due to atomic (H I) and molecular hydrogen 
(H2) along the sight-line often complicates the far-UV spec- 
tra of these objects. Toward a CSPN, this sight-line material 
typica lly consists of interstellar and c ircumstellar components 
(e.g.. lHerald & Bianchill2002l 12004a b). Material comprising 
the circumstellar H I and H2 presumably was ejected from the 
star earlier in its history, and is thus important from an evolu- 
tionary perspective. 

The H2 transitions in the FUSE range are from the Lyman 
(B l E + u -X l Ep and Werner (C 1 n±-Z 1 S+) series. We (HB04) 
found that, for our previous LMC CSPN, a significant amount 



of hot (i.e., T > 2000 K) H2 lies along the sight-lines, presum- 
ably associated with the nebula. At such temperatures, higher 
vibrational states become populated and the absorption pat- 
tern is much more complex, obscuring the entire far-UV con- 
tinuum spectrum (see Fig. 5 of HB04). With care, the effects 
of such a gas can be accounted for and strong stellar features 
can be discerned, but weaker stellar features will be unrecov- 
erable. 

We have modeled the H2 toward our sample CSPN in the 
following manner. For a given column density (AO and gas 
temperature (T), the absorption profile of each line is calcu- 
lated by multiplying the line core optical depth (to) by the 
Voigt profile [H(a,x)] (normalized to unity) where x is the 
frequency in Doppler units and a is the ratio of the line damp- 
ing constant to the Doppler width (the "b" parameter). The 
observed flux is then F„b s = exp[-T(,H(a,x)] x F mtrms i C . We 
first assume the presence of an interstellar component with 
T = 80 K (correspond ing to the mean temperature of the ISM 
— iHughes et al.l 19711) and Vf M rfe=10kms 1 . The column den- 
sity of this interstellar component is estimated by fitting its 
strongest transitions. If additional absorption features due 
to higher-energy H2 transitions are observed, a second (hot- 
ter) component is added, presumed to be associated with the 
nebula. We thus velocity-shift the H2 absorption features 
of the circumstellar component to the radial velocity (from 
Table [3]) of the particular star. HB04 found that, typically, 
the absorption pattern of the circumstellar hot H2 was ad- 
equately fit by absorption models with T e ff ~ 2000 K and 
logN ~ 16.7 cm" 2 , and the reader is referred to that work for 
a more detailed discussion of the H2 fitting. We note that 
our terminology of "circumstellar" and "interstellar" compo- 
nents is a simplification, indicating that we assume the "cool" 
component to be mostly interstellar and the "hot" component 
mostly circumstellar. However, the column density derived 
for the cooler "interstellar" component may also include cir- 
cumstellar H2. Our derived parameters for sight-line Hydro- 
gen are shown in Table[5] Note for some objects (e.g., LMC- 
SMP 29 and SMC-SMP 1), we have not attempted to fit the 
H2 spectrum, due to the lack of useful FUSE data. If a number 
appears in the table without uncertainties, that value has been 
assumed rather than derived. Molecular hydrogen disassoci- 
ates around T ~ 2500, so such high temperatures are signs of 
non-equilibrium conditions (note our models assume thermal 
populations). Such H2 characteristics are observed in shocked 
regions. A shocked region is likely complex, containing areas 
of H2 of differing characteristics. Because of the high com- 
plexity of the absorption pattern at these hot temperatures, we 
did not attempt to fit every feature, but the absorption pattern 
in the far-UV as a whole, to a point where the stellar features 
could be identified and modeled. 

When possible, H I column densities are determined from 
the Lya and Ly/3 features, which encompass both the inter- 
stellar and possible nebular components. The absorption pro- 
files of H I are calculated in a similar fashion as described 
above for H2. 

We also list in Table [5] the reddenings implied by our 
measured column densities of H I usi ng the relationshi p 
(W(Hl)/£ B -v) =4.8xl0 21 cm" 2 mag" 1 (iBohlin et al.lll978h . 
which represents typical conditions in the (Galactic) ISM. 

4.2. The Nebular Continuum 

Nebular parameters taken from the literature for the pro- 
gram objects are compiled in Table [3] They include the angu- 
lar size of the nebula 8, the nebular radius r ne b, the expansion 
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velocity v exp , the electron density n e , the electron temperature 
T e , the H/3 flux Fh/3, the Helium to Hydrogen ratio and the 
doubly to singly ionized Helium ratio. Values in bold are our 
derived values (see below). We also list the reddenings deter- 
mined from literature values of the logarithmic extinction at 
H/3 using the relation cnp = 1 .475 £b-v- 

In some cases, the nebular continuum significantly con- 
tributes to the observed spectra at wave lengths > 1200 A and 
must be estimated to fit the data (§ 14. 31 . 

The nebular contin uum emission has be en modeled using 
the code described in Bianch fet al.l (119971) . which accounts 
for two-photon, H and He recombination, and free-free emis- 
sion processes. The computed emissivity coefficient of the 
nebular gas is scaled as an initial approximation to the total 
flux at the Earth by deriving the emitting volume from the 
dereddened absolute H/3 flux. In the case of LMC-SMP 83, 
the nebula was not entirely contained within the aperture used 
for the STIS observations, and so we scaled the continuum 
by an appropriate geometrical factor (4A/ir9 2 where A is the 
area of the aperture in square arcseconds). T e and n e from the 
literature are used as initial inputs, and adjusted if necessary. 
LMC-SMP 1 and LMC-SMP 76 have no published values of 
the He 2+ /He + ratio. For these, we found adopting a ratio of 
to produce satisfactory fits (these are relativity cooler objects). 
The values of Fnp and Cjj/3 are perhaps the most uncertain, and 
we finally scale the nebular flux with the assumption that it is 
responsible for the majority of the flux at longer wavelengths. 
The nebular continuum flux level is further constrained not to 
exceed the P-Cygni troughs in the UV range. Our combined 
stellar and nebular models, along with the observations, are 
shown in Fig. [5] 

4.3. The Central Stars 

To m odel the spectra of the centr a l stars, we used C M- 
FGEN (lHiUier& Milled fl998l [19991 iHillieret alJl2003h . a 
line-blanketed code suitable for an extended, spherically- 
symmetric expanding atmosphe re. We also used the h ydro- 
static structures from TLUSTY dHubenv & Lanzll 19951) mod- 
els as input to the CMFGEN models. TLUSTY calculates 
the atmospheric structure assuming a plane-parallel geometry. 
The detailed workings of the code are explained in the refer- 
ences above, here we give a brief descriptions of the salient 
features. 

For CMFGEN models, the fundamental photospheric/wind 
parameters include T e ff, R*, the mass loss rate M, the ele- 
mental abundances, the velocity law and the wind terminal ve- 
locity, Vqo. is taken to be the inner boundary of the model 
atmosphere (corresponding to a Rosseland optical depth of 
~ 20). The stellar temperature 7 1 * is related to the luminosity 
and radius by L = AnR^a r* 4 , and the effective temperature 
( T e ff) is similarly defined but at a radius corresponding to a 
Rosseland optical depth of 2/3 (/?2/3)- The luminosity is con- 
served at all depths, so L = A-nR^j^a T e jj 4 . 

We assumed what is essentially a standard velocity law for 
the stellar wind v(r) = v oc (l -r^/r) 13 where ro is roughly equal 
to R*. The choice of velocity law mainly affects the profile 
shape, not the total optical depth of the wind lines, and does 
not greatly influence the derived stellar parameters. Once a 
velocity law is specified, the density structure of the wind 
p(r) is then parameterized by the mass-loss rate M through 
the equation of continuity: M = An R* 2 p(r)v(r). However, 
one actually can only derive M S m=(M c [l ^ff) from the mod- 
els, where M sm and M c ; are the smooth and clumped mass- 



loss rates, and the degree of d umpiness is parametrized by the 
filling factor (/). Both theory (lOwocki etal.ll 19881 1 1994 and 
observations indicate that radiation driven winds are clumped 
to some extent. Unless otherwise noted, M refers to M sm 
throughout this paper. 

It has been found that wind models with the same tnms- 
formed radius R, [oc R*( Voo/M) 2/3 ] (iSchmutz et all 1 19891) 
and Vqo have the same ionization structure, temperature strat- 
ifi cation and spectra (as i de from appropriate scalings with R* 
— ISchmutz et al.lll989l:lHamann et al.lll993l) . Thus, once the 
velocity law and abundances are set, one parameter may be 
fixed (say /?») and parameter space can then be explored by 
varying only the other two parameters {e.g., M and T e ff). For 
opacities which are proportional to the square of the density, 
the optical depth of the wind scales as oc 7?,~ 3 , so R, can be 
thought of as an optical depth parameter. 

In the CMFGEN models presented by HB04, a constant 
scale height (h) was adopted, which connects the spherically 
extended hydrostatic outer layers to the wind. The gravity is 
related to the scale height as h oc g~ l . Our models here differ 
in this respect, as instead of adopting a constant scale height, 
we have used the hydrostatic structures as generated from our 
TLUSTY m odels in the C MFGEN models, in the manner as 
described by iBouret et alJ (120031) . Because of the severe H2 
absorption in the far-UV, and the masking of the stellar flux 
by the nebular continuum at longer wavelengths, the gravity 
cannot be well constrained in our analysis. Thus, we adopt 
gravities typical of CS PNe of the given T f ff, gu ided by the 
\ogg- T e ff relations of IVassiliadis & Wo od ( 1994) (hereafter, 
VW94) (e.g., \ogg = A.l for T eff ~ 55 kK, logg = 6.0 for 
T e ff ~ 100 kK, etc.). The wind features (used as diagnostics) 
are not very sensitive to the gravity, and typically the observed 
wind features can be fit with models of the same T e ff but with 
gravities differing by over a magnitude. Line profiles in the 
optical, which would provide good diagnostics of the grav- 
ity, are not available for our sample (the G750L STIS data 
of LMC-SMP 83 does cover the optical range, but due to its 
emission line spectrum, provides no useful gravity diagnos- 
tics). 

4.3. 1. Abundances and Model Ions 

We note here that throughout this work, the nomenclature 
Xj represents the mass fraction of element i, and 'X Q " de- 
notes the solar value. Following HB04, we adopted two model 
abundance patterns for the MC CSPN. The first, appropriate 
for "H-rich" CSPN, have solar abundances for H and He but 
assume LMC/SMC metallicities of z = 0.4/0.1 Z for the 
metals. For the second, corresponding to "H-deficient" (or 
H-poor) CSPN, we have assumed no hydrogen and an abun- 
dance pattern of X H JX C IX N IX = 0.55/0.37/0.01/0.08 with 
LMC/SMC metallicities of z = 0.4/0.1 Z for the other el- 
ements (e.g., S, Si, & Fe). The abundances for the H-rich 
and H-deficient grids are summarized in Table [4] H-deficient 
abundances (which are enriched in oxygen and carbon) are 
signs of He-burning. Normal, "H-rich" abundances can be 
seen when the star is either H or He-burning. 

We mention here that there are no diagnostics in the far-UV 
or UV spectra of our objects that allow us to make a state- 
ment regarding their hydrogen/helium ratio. The handful of 
wind lines upon which we are basing our parameter determi- 
nations do not allow us to make firm statements regarding the 
specific abundances of a given element, but in most cases we 
are able to determine which grid is more appropriate based on 
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the strength of carbon and oxygen features. 

For the model ions, both TLUSTY and CMFGEN utilize 
the concept of "superlevels", whereby levels of similar ener- 
gies are grouped tog ether and treate d as a single level in the 
rate equations (after I Andersonl[l989l) . Ions and the number of 
levels and superlevels included in the model calculations, as 
well references to the atomic data, are given in the Appendix 

(8EJ. 

4.3.2. Diagnostics 

Our sample covers a wide range of temperatures and mass- 
loss rates. As different diagnostics are used in different pa- 
rameter regimes, we shall discuss the diagnostics used in a 
case-by-case basis (§ [5]l and here we discuss only the more 
general diagnostics. 

For stars showing strong wind features in their spectra, the 
terminal wind velocity ( v^) can be estimated from the blue 
edge of the P-Cygni absorption features. Strong P-Cygni pro- 
files unobscuredby H2 absorption or nebular emission are rare 
in the far-UV spectra of our sample. We mostly used C IV 
AA1548-51 (when present) for our initial estimate of Voo, and 
adjusted this value to match the wind lines. Due to lack of 
suitable diagnostics, we did not attempt to rigorously con- 
strain the degree of clumping in the winds of our sample. In 
the parameter regime of our objects, the model spectra are 
mainly insensitive to /, except for O VI AA1032,38 in the hot- 
ter models (T eff > 50 kK) and P V AA1 1 18,28 for cooler mod- 
els (HB04). We note here that the smooth mass-loss rate M sm 
is an upper limit of the actual mass-loss rate. 

For the other parameters, we adopted the following method. 
We compared the observed spectra with the theoretical spec- 
tra of our model grids to determine the stellar parameters T e ff 
and R, (in the case of CMFGEN models) or T e jj and logg 
(in the cast of TLUSTY models). Far-UV and UV flux levels 
(corrected for the nebular contribution) were used to set 
(since the distance is known), and thus M and L can be de- 
rived from R, . Uncertainties in the quoted parameters reflect 
the range of acceptable model fits. 

Typically, one determines T e ff by the ionization of the 
CNO elements. For these objects, we typically use C III/C IV 
and/or O IV/O V/O VI to constrain T e ff- It is desirable to use 
diagnostics from many elements to ensure consistency. The 
presence or absence of some features offer further constraints. 
For example, for mass-loss rates in the regime of our stars, 
Si IV AA1394,1402 disappears at temperatures T eff > 50 kK, 
S IV AA1062, 1075 at ~45 kK, and P V A A 1 1 18,28 at -60 kK. 
We generally use all wind lines to constrain M. 

5. STELLAR MODELING RESULTS 

Our derived central star parameters are summarized in Ta- 
ble [6] We also indicate the model abundances used for each 
object, and list the reddenings determined by fitting the multi- 
ple components to the far-UV/UV spectra. We now discuss 
the modeling diagnostics and results for the individual ob- 
jects. 

5.1. Results for LMC CSPN 
LMC-SMP 76 (Fig.® 

The far-UV and UV spectra of this object are rich in wind 
features: S VI AA933,44, C III A977 & 1 175, P V AA1 1 18,28, 
O VI AA1032,38, N V AA1238,43, O V A1371, and C IV 
AA1548,51. The O VI doublet has particularly strong ab- 
sorption (but some of the emission seen includes O I airglow 
lines). 



The strong carbon and oxygen spectrum exhibited by this 
object are best fit by models from our H-deficient grid (im- 
plying He-burning). In the model shown ( T e ff = 50 kK, 
M = 3.15 xl0~ 8 M Q yr _1 ) most of the features mentioned 
above are fit adequately, except P V AA1 118,28 (which is 
weak) and C III Al 175 (which is strong). Adjusting the 
temperature/mass-loss rate within the quoted uncertainties al- 
leviate these problems, at the expense of degrading the fits 
of some of the other diagnostics. In cooler models the C III 
features becomes too strong, unobserved Si IV AA1394,1402 
appears, and N V AA 1238-43 disappears. Hotter models 
( T eff > 60 kK) lose the C III lines. Models fail to duplicate 
the strong O VI line in any case, which requires T e ff > 70 kK 
before it begins to be fit adequately. This problem is allevi- 
ated by including X-rays in the model atmosphere (discussed 
below). Based on the C III and C IV wind lines, we derive a 
terminal wind velocity of Voo = 1250 ±250 kms" 1 . 

Our parameters for this star fall between the (M,-,M/) = 
(1.0, 0.578) M Q and (1.5, 0.626) (He-burning) tracks (z = 
0.008) of VW94. This, combined with a better fit of the ob- 
servation with models from our H-deficient (i.e., C and O- 
enriched) grid, cause us to conclu de LMC-SMP 76 t o be a 
He-burner. This is in contrast to Do pita et~aT1 (1 19971) . who 
concluded this object to be H-burning from comparison with 
evolutionary tracks based on their higher derived luminosity 
of logL = 3.735 Lq. They also inferred a very high nebu- 
lar carbon abundance for this object, which is c onsistent with 
our fin dings that the wind is carbon-enriched. Bianc hTet al.l 
d 1997b determined T eff = 52 kK, v^ = 1500 kms" 1 for this 
object by modeling its UV and optical spectra, basically con- 
firmed by our analysis. 

The FUSE spectra of this star show the signature of hot 
molecular hydrogen, which severely attenuates the stellar 
flux. We have fit the absorption pattern with circumstellar 
H2 component having a temperature of T ~ 3000 K. 

As mentioned above, our models failed to fit the O VI 
AA1032,38 line simultaneously with the other, lower- 
ionization diagnostics. It is known that X-rays, arising 
from shocks believed to develop from instabilities in the 
radiationally-driven wind, are needed to explain observed 
O VI l ine profiles in O and ear l y B-stars dBianchi & Garcia] 
12001 iGarcia & Bianchll 120041: iBianchi et al.l l2006h . We 
thus calculated some models which include X-rays in 
the wind. CMFGE N implements X-rays as described in 
iMartins et al.1 (120051) . In summary, shock sources are as- 
sumed to be distributed throughout the atmosphere once the 
wind achieves a certain ve locity, with emissivities taken from 
iRavmond & Smith! d!977l) . The shock parameters include the 
shock temperature (which controls the energy distribution of 
the shocks), the volume filling factor (which sets the level 
of emission), and the velocity of the wind where the shocks 
are switched on. For the starting velocity, we have used 
— 400 kms" 1 , as radiative shocks are believed to become im- 
portant when the wind achieves this velocity. We adopted a 
shock temperature o f 300 kK, typical of high energy pho - 
tons in O type stars (ISchulz et al.ll2003[ ICohen et al.l 120031) . 
and adjusted the volume filling factor until adequate fits of 
the O VI feature were achieved. CMFGEN outputs the to- 
tal X-ray luminosity, and the X-ray luminosities in the en- 
ergy range above 1 and 0. 1 ke V, for the fluxes both absorbed 
by the wind, the X-ray flux which escapes. For our adopted 
parameters, the X-ray luminosity above 0.1 keV is about 5 
times that above 1 keV. We found that for log L x / Lboi > -9.0, 
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the strength of the absorption and the red emission compo- 
nents of the O VI P-Cygni profile were adequately duplicated 
(see example Fig. |7J, although the blue emission compo- 
nent is too strong (most probably due to our simple velocity 
law). For greater L x , the emission feature was a bit stronger 
than observed. Other spectral features did not change until 
log L x / Lboi — -6.0, at which point N IV A955 decreased in 
strength. The O VI feature did not change with respect to 
our model without X-rays for log L x / Lboi < -10.0. Note that 
these X-ray luminosities represent the escaped X-ray flux, i.e., 
what would be observed. The actual X-ray energy in put into 
the wind is ~ 500 times more. iGuerrero et al.1 (1200 ll) discuss 
the possibility that shocks in the stellar wind may explain the 
observed, point-source X-ray fluxes (log L x / Lboi < -7.0) of 
the central star of the Cat's Eye nebula (NGC 6453). 
LMC-SMP 1 (FigM 

The spectrum of LMC-SMP 1 is similar to that of LMC- 
SMP 76, but with the wind features having a lower terminal 
velocity (Vqo < 800 kms -1 ) and weaker oxygen signature. 
We can achieve adequate fits using either models from the 
H-rich or H-poor grid, with the former requiring higher mass- 
loss rates. We list both models in Tabled but only show the 
H-poor model in the figures. The higher temperature limit 
( T eff = 70 kK) fits the O VI AA1032,38 absorption trough bet- 
ter. For the lower temperatures ( T e ff = 55 kK) and mass-loss 
rates, the C III A 1175 feature is too strong compared with 
the observations. In the UV, the N V AA1238,43 and C IV 
AA1548,51 are adequately fit with models in the quoted pa- 
rameter range, but the O V A 1371 feature is always too strong 
(HB04 had similar problems with this line with some of their 
sample) in the case of either abundance pattern. The O IV 
AA1339,43 feature is similarly too strong in the H-poor model 
spectra, while the carbon features are generally better-fit than 
with our H-rich model grid (probably indicating our adopted 
oxygen abundance is too high). The FUSE spectrum shows 
the absorption pattern from hot (T ~ 3000 K) molecular hy- 
drogen, which severely depresses the stellar flux. In compar- 
ison with LMC-SMP 76, this object has a higher luminosity 
and effective temperature, but a lower Vqo (mass-loss rates 
are comparable). This could be explained by our finding that 
the wind of LMC-SMP 76 seems more obviously chemically 
enriched, and thus a higher opacity leading to a more efficient 
transfer of radiation-to- wind momentum. 
LMC-SMP 83 (FigM: 

LMC-SMP 83 is a bit of an oddball, for which the classifi- 
cation as a CSPN is a matter of current debate. The object 
is notable for multiple reasons. Between 1993 and 1994, its 
luminosity increased from logL = 4.6 to 5.6 (L ). The qui- 
escent luminosity is already extre me for a CSPN (most LMC 
CSPN have logL < 3.9 L Q , e.g., iDopita & Meatheringhaml 
Il991a[b|), but is less luminous tha n any known massive WN 
(e.g.. lHamann & Koesterkelll998l) . Furthermore, it is unique 
in that it displays a [WN] spectrum (other LMC CSPN which 
have WR-type spectra show [WC] spectra). It has been the 
subject of a long-term monito ring campaign and UV spectra 
exist spanning many epochs. Ham ann et al.l J2003) (hereafter, 
HPGR) performed detailed UV and optical spectral model- 
ing across these epochs, and found the following parameters 
to be constant: T* = 112 ±20 kK, X H = 0.2, X N = 3 xlO -3 , 
Xc < 1 xl0~ 4 , Xf e = 2 xlO -4 . The increase in luminosity dur- 
ing outburst was apparently caused by an increase in mass- 
loss rate, which increased the effective radius. The nature of 
LMC-SMP 83 is debatable, and HPGR discuss many alter- 
natives to the CSPN classification (e.g., single massive star, 



low-mass binary system with mass transfer, double degener- 
ate merger). The main argument against the CSPN classifica- 
tion is its surface nitrogen-rich abundances, which are difficult 
to explain in the CSPN context. 

In Fig. [10] we show the absorption pattern of the interstel- 
lar hydrogen absorption component, and an additional hot cir- 
cumstellar component applied to a flat continuum. There does 
appear to be an overall signature of hot H2 in the data, al- 
though many minor features are not well fit. The H2 absorp- 
tions are significant, and no obvious strong stellar features are 
apparent in the FUSE observations. Because the FUSE range 
offers no further stellar diagnostics, we have not performed 
much further modeling of this object beyond that of HPGR. 

The FUSE flux levels are a factor of about 1 .5 higher com- 
pared to the STIS data. LMC-SMP 83 was observed by STIS 
in 1999 and again in 2000, and comparison of these two data 
sets shows the 1999 continuum level in the 1250-1550 A re- 
gion to be about 95% that of the 2000 level measured 13 
months later (the FUSE data was gathered abo ut 1.5 yr af- 
ter the 2000 STIS observations). As noted by lHamann et alj 
d2003l) . there is very little change in the actual spectral fea- 
tures within this time. For this work we have used the 1999 
STIS data set, because its signal to noise is superior to that of 
the 2000 data set, and scaled the FUSE data to the STIS flux 
levels. 

We first computed a model with parameters similar to 
what HPGR found for the most recent epoch (2000): Voo = 
1600 kms" 1 , R*=0.5 R Q , T eff = 95 kK, log M = -5.67 M , 
and the constant parameters quoted above (abundances and 
I 1 ,). However, rather than using their reddening (a Galactic 
foreground reddening of Lb-v = 0.03 mag and an LMC red- 
dening of Lb-v = 0.13 mag), we used a smaller amount of red- 
dening, to produce better agreement over our expanded wave- 
length range. The use of this less severe reddening results 
in a slightly smaller derived radius (and corresponding mass- 
loss rate) than those of HPGR. This model produced several 
He II features in the FUV which do not appear to be present in 
the observations. However, it is possible that severe absorp- 
tion from a hot circumstellar H2 gas (§ 14. Il l could attenuate 
the He II features to agree with the observations, as shown in 
Fig. [9] Another possibility is that the mass-loss rate has since 
declined to the point these features no longer appear in emis- 
sion. Note that because we have essentially adopted the stellar 
parameters of this object from HPGR and not performed any 
further stellar modeling, we have omitted uncertainties for its 
stellar parameters in Table [6] 
LMC-SMP 50 (Fig. QTJ 

The photoionization temperature for this star is 100 kK (see 
Table|7]l, and we find that a H-rich model with T e ff— 1 10 kK, 
logg ~ 6.0, logM = -8.2 M , ^ 500 kms" 1 , R* ~ 
0.21 R does an adequate job at fitting the most apparent 
wind feature displayed — that of the O VI doublet (see Fig.[TJ. 
These parameters lie slightly below the M,„„ =2.0 M , M = 
0.67 M track of VW94. The rest of the theoretical spectrum 
is relatively featureless, however, whereas the FUSE and FOS 
display some complexity. Many of the absorption features in 
the shorter wavelength region of the FUSE spectrum can be 
attributed to molecular hydrogen, as shown in the figure. The 
poor quality of the observations make firm identification of 
stellar and interstellar features difficult. For example, the FOS 
spectra show a P-Cygni-like feature at A 1335, which might be 
a C II feature of stellar origin (indicating a much cooler object, 
possibly a companion) or possibly a circumstellar absorption 
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feature in a relatively noisy spectrum that mimics an actual P- 
Cygni line. The lack of strong, clean stellar features prohibit 
a more rigorous determination of stellar parameters. 
LMC-SMP 29 

Along with SMC-SMP 22, the FOS spectra show show 
high-ionization lines (e.g., N V AA1238,43) and have extreme 
photoionization temperatures (T ~ 200 kK). Unfortunately, 
problems during the FUSE observation (§|2]i resulted in unre- 
liable FUSE data for this object, which we have not attempted 
to model. We have also made no attempt to fit the interstellar 
Lya absorption feature. 

5.2. Results for SMC CSPN 

SMC-SMP 3 (Fig. [/2j 

There are no obvious wind features in the FOS spectrum 
of this object, and the FUSE spectrum only obviously shows 
the O VI resonance line. The FUSE spectrum is also contami- 
nated by several strong airglow features, especially at shorter 
(< 1000 A) wavelengths, and it shows a moderate IT signa- 
ture, which we have fitted. This star has a relatively high pho- 
toionization temperature of T e ff = 92 kK (Table |7). Given 
the paucity of wind diagnostics, we can only speak in general 
terms about its parameters. Models with 90 < T e ff < 100 kK 
with -8.0 < log M < -7.0 M Q and v^ > 2000 kms" 1 ade- 
quately reproduce the O VI feature. Cooler temperatures and 
higher mass-loss rates typically show unobserved spectral fea- 
tures, and lower mass-loss rates result in a O VI weak com- 
pared to observations. Lower Vqo results in both components 
of the O VI doublet being resolved — contrary to observa- 
tions. However, the P-Cygni edge of the O VI feature is ob- 
scured by absorption due to Ly/3, so we cannot make a more 
precise measurement of Vqo. We conclude by saying a H-rich 
model with temperatures consistent with the photo-ionization 
temperature with M and Voo within the previously stated lim- 
its fit the FUSE and FOS spectra adequately, with the caveat 
that these parameters are constrained mainly by the O VI dou- 
blet and the absence of other strong features. For comparison 
with evolutionary models, a T e ff = 95 kK model scaled to the 
observed FOS flux levels would have ~ 0.32 R , logL ~ 
3.8 L Q . Comparison with the Z = 0.004, M Mt = 1.5 M Q , 
M = 0.640 M© track of VW94 implies a post-AGB age of a 
few kyr. 

For this object, as well as LMC-SMP 50, we have relied 
mainly on the O VI doublet as a diagnostic. In our discus- 
sion of LMC-SMP 76, we showed this line could be heavily 
influenced by X-rays, which might add additional uncertainty 
to parameters derived using this line. However, photoioniza- 
tion models indicate that SMC-SMP 3 (and LMC-SMP 50) 
are considerably hotter than LMC-SMP 76 (T eff > 90 vs. 
~ 50 kK). The absence of low-ionization spectral features 
supports the idea that these two objects lie in this high tem- 
perature regime. Preliminary models computed by us indicate 
that, in this regime, the high temperatures alone are adequate 
to reproduce the observed O VI profile, and including reason- 
able amounts of X-rays in the model atmospheres does not 
alter this profile significantly. We plan a more detailed inves- 
tigation of the role of X-rays in the winds of CSPN to be the 
subject of a future paper. 
SMC-SMP 1 (Fig.\ij) 

Due to uncertainties with the FUSE observations, we have 
determined the characteristics of this object by mainly re- 
lying on wind features in its UV spectrum, using Si IV 
AA1394,1402, and C IV AA1548,51 as diagnostics. We have 



also used the far-UV C III Al 175 feature for guidance. These 
features indicate a relatively low terminal wind velocity of 
Vqo — 500 kms" 1 . Because of the limited diagnostics, we 
have assumed normal, "H-rich" SMC abundances. Models 
of T e ff ~ 37 kK produce acceptable fits to the UV spectrum 
(cooler temperatures weaken C IV AA1548,51 unacceptably, 
and hotter temperatures prohibit the Si IV AA1394,1402 and 
C III Al 175 features). However, we are unable to attain ac- 
ceptable fits of the FUSE spectrum with our stellar models, 
even when incorporating various hydrogen absorption mod- 
els. The limited wavelength coverage prevent a determination 
of £b-v from the spectral fitting, s o we have simply adopted 
the va lue of £b-v = 0.16 (from Meatheringham & Dopita 
Il991al see Table [3). We note that in Fig. [T3l what appears 
to be strong stellar He II emission feature at A 1640 is actu- 
ally multiple narrow emission features (one of which being 
nebular He II) blended together as a result of the binning for 
clarity. 

SMC-SMP '22 (Fig.\]J§ 

This star has a very high photoionization temperature of 
T e ff = 200 kK (Table |7). However, we are able to fit the 
FUSE spectrum well with a SMC-metallicity TLUSTY model 
of T e ff = 25 kK (logg = 3.7). Most of the detailed struc- 
ture is from Fe III, which is present in the atmosphere for 
T eff < 30 kK With z = 0.1 Z Q , the C III A1175 feature is 
a bit weak. Using LMC metallicity (z = 0.4 Z Q ) the line is 
better fit, but a higher temperature ( T e ff = 25) is required to 
match the Fe III spectrum. We are able to fit its FUSE spec- 
trum without invoking a hot circumstellar molecular gas as 
found in other MC CSPN. 

The photometric SMC catalog of Zari tskv et ail d2002l) lists 
an object 17.4 arcseconds away from our target (coordinates: 
58 36.77 -71 36 6.23) with B-V = -0.26, U-B = -0.77. These 
colors are consistent with those of an early B-type star, which 
would have T e ff ~ 20 kK. A main sequence star of this sub- 
type would have R ~ 4 R© and logg ~ 4. Using a TLUSTY 
model with these parameters, and scaling the flux to the dis- 
tance of the SMC (assuming a foreground reddening of 0.05) 
results in a far-UV flux level of ~ 1 xlO" 14 erg s" 1 cm" 2 " 1 , 
about 80% of the flux levels of our FUSE observations. Scal- 
ing its flux by a factor of ~ 1.5 to match the U-magnitude 
from the catalog (U= 16.38) results in essentially the far-UV 
flux levels of the FUSE observations. Given that the FUSE 
flux levels are several times higher than the FOS fluxes, we 
believe that the most likely explanation is that, due to the sig- 
nificant target drift experienced during the observation of this 
target (§ 12. U . the light from the B-type star entered the FUSE 
aperture to some extent, and the FUSE spectra are dominated 
by this other object, rather than the target. Thus we exclude 
this object from subsequent analysis. 

6. DISCUSSION 

In Table|7] we present a combined summary of our findings 
here with those of HB04. We list temperatures derived fro m 
photo-ionization models (jDopita & Mea theringhaml[r 991alfbl : 
iDopita et al.lfl99 7; Vassil iadis et al.lfl998bl) as well as those 
from our stellar atmospheres analysis. We also list the mea- 
sured terminal velocity of the wind, and the morphology of 
the surrounding PN. In Figs. [TBI and [TBI we plot our sample 
along with that of HB04 with the evolutionary tracks of V W94 
for LMC and SMC metallicities, respectively (both hydrogen 
and helium burning tracks are shown). Most of the the LMC 
sample is seen to lie in the He-burning regime. VW94 found 
similar results when they plotted their tracks with the CSPN 
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samples of Dopita & Meatheringham ( 1991a b) (who derived 
parameters based on photoionization models), and attributed 
this result to the different timescales of the H and He-burning 
tracks. He-burners are expected to take several thousands 
of years to migrate across the interval logL ~ 4.0 — 2.5 L Q , 
whereas the H-burners only take < 1000 years to span the 
same range. The tracks of VW94 apparently indicate that, as a 
consequence of their faster evolution to smaller luminosities, 
the LMC H-burners have already faded below the threshold 
of FUSE, as our sample represents the expected UV-brightest 
CSPN. 

We found evidence of hot molecular hydrogen in the cir- 
cumstellar environments of our previous sample of 7 LMC 
stars (HB04). Here we have made similar findings for our 
targets. This is further indication that hot H2 in the cir- 
cumstellar environment s of PN is not uncommon (e.g., see 
iHerald&Bianchil 120021 l2004bl for Galactic examples). H 2 
may exist in clumps, shielded from the intense UV radia- 
tion fields by neutral and ioni zed hydrogen, as a ppears to be 
the ca se in the Helix nebula (Speck et al. 2002). Spe ck et all 
(2002) suggest that these clumps may form after the onset 
of the PN phase, arising from Rayleigh-Taylor instabilities at 
either the ionizatio n front or the fast wind shock front. The- 
oretical models by Alem an & Gruenwaldl (120041) show that a 
significant amount of H2 can survive within the ionized re- 
gion of PN, and that the ratio of total H2 to hydrogen mass 
inside the ionized region increases with the temperature of 
the central star. A correlation between the presence of H2 line 
emission and the morpholo g y of the PN has bee n noted by 
IZuckerman & Gatleyl (ll988l) : lKastner et al.l (1 19961) (known as 
"Gatley's rule"). 

For the stars analyzed in this paper, our derived effective 
temperatures are mostly in good agreement with those de- 
rived from photoionization models. For SMC-SMP 3, which 
shows sparse spectral diagnostics, a stellar model with a tem- 
perature of about the photoionization temperature (~ 95 kK) 
fits the spectra adequately. The temperature derived through 
spectral modeling (by HPGR) of the enigmatic LMC-SMP 83 
is significantly lower than through photoionization models 
(~ lOOkKvs. 170 kK). 

7. CONCLUSIONS 

HB04 analyzed 6 LMC CSPN, and here we have analyzed 
FUSE observations of 4 additional LMC and 2 SMC CSPN. 
We determined the stellar parameters using the stellar atmo- 
sphere code CMFGEN to analyze their FUSE far-UV and 
HST UV spectra. In some cases, we also modeled the neb- 
ular continua (if it contaminates the UV spectra) and the 
atomic and molecular hydrogen absorption along the sight- 
line (which severely affects the far-UV flux). 

To summarize, we find: 

• For those objects with photoionization temperatures be- 
tween 30 < T < 100 kK, the temperatures determined 
via stellar modeling are in rough agreement. For LMC- 
SMP 83, our derived effective temperature of T e ff ~ 
95 kK, is significantly below the photoionization tem- 
perature of T = 170 kK. However, this unique object 
has shown substantial changes in its spectra (and, con- 
sequently, its parameters) over the past two decades. 

• Most objects often have absorption features in the far- 



UV which could be attributed to very hot (T > 2000 K) 
H2 in their circumstellar environment. These temper- 
atures may be due to the proxim ity of the nebulae to 
the st ar (see theoretical work of lAleman & Gruenwaldl 
2004)). or perh aps due to shocks (as suggested by 
Speck et alJl2002l) . 

• The majority of our objects have Voo < 1000 kms" 1 . 
Two of the exceptions show signs of chemically en- 
riched winds (LMC-SMP 83 & LMC-SMP 76). For 
the other (SMC-SMP 3) the poor data quality prohibit 
a definitive statement about wind abundances from this 
spectral analysis. 

• We find indications of X-rays in the winds of LMC- 
SMP 76. We are able to attain good fits consis- 
tently for the stronger wind features, except for the 
far-UV O VI AA1032,38 doublet, which is too weak 
in our models unless X-rays are included in the cal- 
culations. Although it has been shown that X-rays 
are n eeded to replicate this feature for O and early Bi - 
sters (IBianchi & Garcia 2002t iGarcia & Bianchill2004fc 
iBianchi et alj 120061) . this is the first case of a similar 
phenomena in CSPN winds. 

• These FUSE observations add two cases (LMC- 
SMP 83 and SMC-SMP 22) to the small number of 
CSPN displaying nebular O VI AA1032,38 (the first 
case, LMC-SMP 62, was reported by HB04). 

By virtue of their membership of the Magellanic Clouds, 
the uncertainties in their distances is small, which carries over 
to the derived stellar parameters. This is a great advantage 
over Galactic CSPN, where the distance is the largest source 
of uncertainty in the analysis. 

Our FUSE observations have allowed us to derive a set 
of stellar and wind parameters for CSPN of the Magellanic 
Clouds that are unhampered by the distance uncertainties that 
plague equivalent Galactic studies. In addition to revealing 
the flux of the hot star (which is obscured at longer wave- 
lengths by the nebular flux), the far-UV observations also in- 
dicate that many of these young CSPN have hot molecular hy- 
drogen in their nebulae, and are thus also giving insight into 
their circumstellar environments and histories. 
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APPENDIX 
APPENDIX: MODEL ATOMS 

For our TLUSTY models, the atomic data used come from TOPBASE, the data-base of the Opacity Project dCunto et alj 
1993). Hydrogen and Helium were treated in NLTE, and the other included species (C, N, O, Ne, Na, Mg, Si, S, Ar, Ca, and Fe) 
were allowed to contribute to the total number of particles and charge but their opacity contribution was neglected in the model 
atmosphere calculation. 

Ions and the number of levels and superlevels included in t he CMFGEN model calculations are listed in ab le [8] The atomic 
data come from a va r iety of sources, with t he Opacity Proj ect dSeatonlll987l : IOpacitv ProiectT eam 19951 1 19971) . the Iron Project 
dPradhan et alJI 19961: iHummer et aill 19931) . iKuruczl (1995b 1 and the Atomic Spectra Database at NIST Physical La boratory be- 
ing the principal sources. Much of the Kurucz atomic data were obtain e d directly from CfA (jKurucz 1988, 2002). Ind ividua l 
sources of a t omic d ata include the fo ll owing: IZhang & Pradhanl d 19971). iBautista & Pradhanl d!997l) . iBecker & Butlerl d 19951) . 
iButler et alj d!993l). iFuhr et all d!988l). iLuo & Pradhanl d!989l). iLuo et alj d!989l). Me ndoza (1983. 1995, private communica- 
tion). iMendoza et all d 19951). iN ussbau mer & Storey! d!983l 11984 . iPeach et all d 19881) . Storey (1988, private communication), 
iTully et ail (Il990h . and lWiese et alJ (1 1 966l IT9 69 ) . Unpublished data taken from the Opacity Project include: Fe VI data (Butler, 
K.) Fe VIII data (Saraph and Storey) and C. Mendoza (Fe IX, Fe X). 
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FIG. 1 . — FUSE spectra of our LMC sample, velocity shifted to the rest frame of the central star (Table|3j, and re-binned to a resolution of 0.10A. The more 
prominent stellar features and nebular emission lines are marked by black labels, airglow features are marked by the blue labels. Most objects seem to display 
wind lines to some extent. Nebular O VI emission appears in the spectra of LMC-SMP 83. 
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FIG. 2. — FOS and STIS (for LMC-SMP 83) spectra of our LMC sample, velocity shifted to the rest frame of the central star (Table|5J, and re-binned to a 
resolution of 0.25A. The more prominent stellar features and nebular emission lines are marked with black labels, interstellar absorption features are marked 
with blue/gray labels. Most objects display wind lines to some extent. C IV AA1548-5 1 appears as a P-Cygni profile in most. The spectra of LMC-SMP 50 and 
LMC-SMP 29 show emission lines from nebulae in high states of ionization. 
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FIG. 3. — FUSE spectra of our SMC sample, velocity shifted to the rest frame of the central star (Table[5J, and re-binned to a resolution of 0.10A. The more 
prominent stellar features and nebular emission lines are marked by black labels, airglow features are marked by the blue/gray labels. Note the O VI AA1032-38 
nebular emission features of SMC-SMP 22. 
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FIG. 4. — UV spectra (FOS) of our SMC sample, velocity shifted to the rest frame of the central star (Table[3), and re-binned to a resolution of 0.25A. The more 
prominent stellar features and nebular emission lines are marked with black labels, airglow lines or interstellar absorption features are marked with blue/gray 
labels. 
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FIG. 5. — The MC CSPN spectra (black) are shown along with our stellar (pink/grey) and nebular continuum models (blue dashed). The sum of the models 
is shown by the red line. All models have been reddened with our determined values for £b-v (Table|6j, and the effects of hydrogen absorption (Table|5j have 
been applied. The spectra have been convolved with a 5 A Gaussian for clarity. Note the logarithmic flux scale. 
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FIG. 6. — LMC-SMP 76: The observations are shown (black) along with our stellar model, with and without our hydrogen absorption models applied (pink/gray 
and blue dotted, respectively). The FUSE spectra have been convolved with a 0.6 A Gaussian for clarity. 
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FIG. 7. — LMC-SMP 76: The observations are shown (black) along with two models, one without X-rays (pink/gray) and one with X-rays of strength 
log L x / Lbol ^ —7.0 (blue/dark grey). Without X-rays, the model fails to duplicate the strong P-Cygni absorption trough caused by O VI AA1032.38. Including 
X-rays in the model calculations strengthens the O VI feature while not affecting other diagnostics. 
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FIG. 9.— LMC-SMP 83: Description follows that of Fig.|6] 
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FIG. 10. — LMC-SMP 83: A portion of the FUSE spectrum is shown (black) along with a flat model continuum, both with only an interstellar hydrogen 
absorption component applied (dashed) and with an additional hot circumstellar H2 component applied (green). 
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FIG. 11 .— LMC-SMP 50: Description follows that of Fig.|5] 
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SMC-SMP 1: Description follows that of Fig.|6] 
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FIG. 14. — SMC-SMP 22: Description follows that of Fig. [6] Note: this observation was apparently contaminated by a B-star in the aperture, and so does not 
represent the true spectrum of SMC-SMP 22 (this is the reason for the flux discrepancy between the far-UV and UV spectra). 




FIG. 15. — HR diagram of our LMC sample and that of HB04. Evolutionary tracks (Z=0.008) of VW94 are shown, labeled with their initial and final masses, 
with the H and He-burning tracks denoted by the solid and dotted lines, respectively. Note our adopted model metallicity corresponds to Z=0.006. 
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FIG. 16. — HR diagram of our SMC sample. Evolutionary tracks (Z=0.004) of VW94 are shown, labeled with their initial and final masses, with the H and 
He-burning tracks denoted by the solid and dotted lines, respectively. Note our adopted model metallicity corresponds to Z=0.0016. 



TABLE 1 

FUSE C056 AND D034 PROGRAM STARS, AND FUSE DETECTOR SEGMENTS UTILIZED 



Name 


R.A. 


Dec 


Archive 


Date(s) 






Wavelength Range (A) 






(2000) 


(2000) 


Name 




SiCl 


SiC2 


LiFl 


LiF2 


LMC-SMP 1 


04 38 34.7 


-70 36 43.9 


D0340101 


7/4/03 


920-993 


916-1000 


1000-1081,1095-1136 


1087-1183,1020-1075 


LMC-SMP 83 


05 36 20.80 


-67 18 08.0 


D0340201 


7/6/03 


915-993 


915-1006 


1000-1084,1093-1184 


1020-1075,1088-1183 


LMC-SMP 76 


05 33 56.20 


-67 53 08.4 


D0340301 


4/29/03 


915-993 


915-1013 


1005-1091,1087-1196 


1005-1082,1079-1186 








D0340302 


4/30/03 










LMC-SMP 50 


05 20 51.82 


-67 05 43.1 


D0340401 


7/3/03 


915-993 


916-1006 


1000-1082,1094-1136 


1028-1082,1086-1183 


LMC-SMP 29 


05 08 03.47 


-68 40 16.5 


D0340501 


7/8/03 






1005-1082,1086-1183 




SMC-SMP 1 


00 23 58.80 


-73 38 04.1 


C0560101 


6/1 1/02 






1005-1084,1093-1187 




SMC-SMP 3 


00 34 21.90 


-73 13 20.8 


C0560201 


6/12/02 


915-993 


915-1005 


1105-1184,1093-1155 


1005-1082, 1088-1183 


SMC-SMP 5 


00 41 21.70 


-72 45 18.1 


C0560402 


8/5/02 
















C0560401 


8/3/02 










SMC-SMP 22 


00 58 37.22 


-71 35 48.8 


C0560301 


7/25/02 


915-992 




992-1084,1093-1187 





TABLE 2 

Archive UV data sets (FOS and STIS) 



Name 


Instrument 


Dataset 


Date 


Grating 


A (A) 


LMC-SMP 76 


FOS 


Y1C10B03T 


9/22/93 


G130H 


1182-1600 




FOS 


Y1C10B03T 


9/22/93 


G190H 


1600-2330 


LMC-SMP 83 


STIS 


O5IM01010 


02/04/00 


G140L 


1182-1720 




STIS 


O5IM01020 


02/04/00 


G230L 


1720-3150 




STIS 


O5IM02010 


02/04/00 


G430L 


3150-5200 




STIS 


O5IM02010 


02/04/00 


G430L 


3150-5200 




STIS 


O55R01010 


01/08/99 


G140L 


1182-1720 




STIS 


O55R01020 


01/08/99 


G230L 


1720-3150 




STIS 


O55R01030 


01/08/99 


G430L 


3150-5200 




STIS 


O55R02020 


01/08/99 


G750L 


5265-10260 


LMC-SMP 1 


FOS 


Y3EV0106P 


11/15/96 


G130H 


1182-1600 




FOS 


Y3EV0103P 


11/15/96 


G190H 


1600-2300 




FOS 


Y3EV0104P 


11/15/96 


G270H 


2300-3280 




FOS 


Y3EV0105P 


11/15/96 


G400H 


3280-4822 


LMC-SMP 50 


FOS 


Y3EV0805T 


11/10/96 


G130H 


1182-1600 




FOS 


Y3EV0802T 


11/10/96 


G190H 


1600-2300 




FOS 


Y3EV0803T 


11/10/96 


G270H 


2300-3280 




FOS 


Y3EV0804T 


11/10/96 


G400H 


3280-4822 


LMC-SMP 29 


FOS 


Y2Y00204T 


10/22/95 


G130H 


1182-1600 




FOS 


Y2Y00205T 


10/22/95 


G190H 


1600-2300 




FOS 


Y2Y00206T 


10/22/95 


G270H 


2300-3280 




FOS 


Y2Y00207T 


10/22/95 


G400H 


3280-4822 


SMC-SMP 1 


FOS 


Y1C10103T 


5/3/93 


G130H 


1182-1600 




FOS 


Y1C10104T 


5/3/93 


G190H 


1600-2330 


SMC-SMP 3 


FOS 


Y1C10203T 


11/26/93 


G130H 


1182-1600 




FOS 


Y1C10204T 


11/26/93 


G190H 


1600-2330 


SMC-SMP 5 


FOS 


Y2Y00904T 


6/12/96 


G130H 


1182-1600 




FOS 


Y2Y00905T 


6/12/96 


G190H 


1600-2330 




FOS 


Y2Y00906T 


6/12/96 


G270H 


2300-3280 




FOS 


Y2Y00907T 


6/12/96 


G400H 


3280-4822 


SMC-SMP 22 


FOS 


Y2N30104T 


3/20/95 


G130H 


1182-1600 




FOS 


Y2N30105T 


3/20/95 


G190H 


1600-2330 
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TABLE 3 
Nebular Parameters 



Name 


gaa 
(") 


r bb 

neb 

(pc) 


v„„(OIII/OIir 
(kms 4 ) 


dyn 

(kyr) 


J2 c ([OII/SII]) b 

(10 3 cm 4 ) 


r c (0 +2 /N +1 )* 

(kK) 


iog(F Hp ) h 



£b-v(c///3)* 
(mag) 


He/H r 


He 2+ /He+ 


VJJSR 
(kms 4 ) 


LMC-SMP 1 


0.33 s 


0.04 


17.2/23.1 


2.3 


3.44/3.32 


11.5/- 


-12.46 (-12.5) 


0.11 


0.146 


- e (0.0) 


209.1° 


LMC-SMP 29 


0.68 


0.083" 


35.9/49.3 


2.3 


5.46/4.06 


20.0/- 


-12.71 


0.05 


0.160 


1.53 E 


228.2" 


LMC-SMP 50 


0.64 s 


0.08 


35.0/55.7 


2.2 


5.37/- 


11.2/- 


-12.71 (-12.81) 


0.04 


0.120 


0.23 e 


284.1" 


LMC-SMP 76 


<0.16 


<0.02 


29.0™/32.9 


0.7 


13.7/- 


12.0/- 


-12.54 


0.22 


0.075 


(0.0) 


262.8" 


LMC-SMP 83 


2.34 


0.292" r 


82.9*780.9 


3.4 


2.46/2.22 


17.4/11.7 


-12.65 


0.10 


0.137 


1.15" 


274.5/300* 


SMC-SMP 1 
SMC-SMP 3 
SMC-SMP 5 
SMC-SMP 22 


0.31 
0.49 
0.42 
0.52 


0.045" 
0.072" 
0.061" 
0.075" 


15.4/16.1 

32.9/- 
29.2/43.6 
50.9/51.2 


2.9 
2.1 
2.0 
1.5 


9.59/- 
1.96/- 
3.7/3.34 
5.51/2.33 


10.7/- 
13.5/- 
14.5/- 
26.6/11.6 


-12.77 (-12.89) 
-13.06 (-12.78) 

-12.81 

-12.88 


0.16 
0.03 
0.07 
0.11 


0.080 
0.160 
0.130 
0.135 s 


-70.155 s 
0.0771.66 s 
0.60 E 


m.<y 

98.2' 
\Q\.2' 
144.2/ 



REFERENCES. — (a): Dopit a etal. (1988) (b): Meatheringham & Dopita i 1991a b) (c): Dopita & Meatheringham 1 1991a b) (d): Vassiliadis et al. 1 1998b) (e): Monk et al. 1 1988) (f): 
lDopitaetal]n996l) (g): lHenrv et all 1 19891) (hTlMeatheringham et al. 1 1988) (\) : iDopita et al J fl985 ) (k): Peiia ct all 120041) (s): IShaw et alj 120011) (x): approx of three vals given 
aa Calculated from r neb using LMC/SMC distances of 50.6/60.0 kpc 
bb From their "D(edge)" 

Line fit with two components 
ee From [S II] l Meatheringham & Dopita ll991bT) 
xx Multiple components 
xxx Avg. of 3 components 



TABLE 4 
Abundances of grid models 



Grid Xh Xn e Xc %N Xo %s ' %p %s %Fe 

LMC H-deficient - 0.54 0.37 0.01 0.08 2.80x10^ 2.45x10"* 1.53 xlO" 4 5.44x10^* 

LMCH-rich 0.71 0.28 1.22xl0~ 3 4.40x10^ 3.82xl0" 3 2.80x10^ 2.45xl0" 6 1.53X10" 4 5.44x10^ 

SMC H-rich 0.71 0.28 3.05x10^ 1.10x10^ 9.55 xl0~ 4 7.00x10^ 6.12xl0" 7 3.65 xlO" 5 1.36X10 -4 



TABLE 5 
H I AND H 2 PARAMETERS 





H I' s+C " 




H 2 ' s 




H 7 C "' 


C 




Star 


log/V 


T 


logW 


T 


log/V 


T 


£b-v(H I) 




(cm- 2 ) 


(K) 


(cm- 2 ) 


(K) 


(cm" 2 ) 


(K) 


(mag) 



0.2 
0.1 

0.05 
0.04 



SMC-SMP 1 2L9 80 ' - ' — 

SMC-SMP 3 21.7 80 19.40+° 3 2.1 80 16.70+Jj' 3 2000±1000 1.0 

SMC-SMP 22 20.00+° 3 80 16.70^ 4 80 - - 0.02 



LMC-SMP 76 21.00; 

LMC-SMP 1 20.70; 

LMC-SMP 50 20.40; 

LMC-SMP 83 20.80 
LMC-SMP 29 



80 18.70+0- 3 80 

80 18.00 ±1.0 80 

80 19.7+°j 80 

80 16.70^ 80 



16.70+°° 2000 ±1000 

17.00™ ; * 3000 ±1500 

16.70+21 3000 ±1500 

16.70+° :3 2000 ±1000 



TABLE 6 
Derived Stellar Parameters 



Name 


Model 
Abundance 


Teff 

(kK) 


logL 

(L ) 


(Ro) 


log M 

(M Q yi-') 


Voo 

(kms" 1 ) 


(cgs) 


£ B -v 
(mag) 


LMC-SMP 83 
LMC-SMP 1 
LMC-SMP 1 
LMC-SMP 76 
LMC-SMP 50 
LMC-SMP 29 


LMC 83 
H-poor 
H-rich 
H-poor 
H-rich 


95 
65+ 5 
65+ 5 
50±5 
110± 10 


4 49+o.iv 
3 61^™ 

T. 46+0-08 

3.80$* 


0.47 ±0.1 
70+° 10 
66+ 8 '0 

0.67 ±0.1 
22+ 03 


-5.71+ 013 
—7 52+ '- 2 ^ 
_ 7 Af.iftit, 

IMa -0.2l 
_7 4Q+0.37 

o 7+0.30 
°-^-0.20 


1600 
800 ±100 
800 ±100 
1250 ±250 
500 ±100 


6.0 
5.0 
5.0 
4.7 
6.0 


0.08 
0.13 
0.13 
0.03 
0.04 
0.05 


SMC-SMP 1 
SMC-SMP 3 


H-rich 
H-rich 


37+ 3 
95±5 


3 5+ 

3.S_ 02 


1.3 ±0.2 
0.32 ±0.05 


-7.31±0.10 
-7.8 ±0.1 


500 ± 100 
> 2000 


4.9 
5.7 


0.16 
0.03 
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TABLE 7 
Summary Table 



Name 


T aa 


phot 


v "" 


Morphology 




(kK) 


(kK) 


(kmr 1 ) 




LML-SMH 2 


38 


39 


700 


spnencal , Round 


t A-rr^ ca/tp a< 

LlVlL^-iJVljr oj 




4j 




Kouna 






127 


1 000 


Bipolar w/central condensation'"* 


LMC-SMP 76 


50 


50 


1250 


Spherical'" , Round 




5< 


46 


1 000 


Bipolar w/central condensation'' 


LMC-SMP 23 


60 


65 


1100 


Ring 6 , Round" 


LMC-SMP 1 


65 


66 


800 




LMC-SMP 61 


70 


59 


1300 


Spherical* 


LMC-SMP 50 


110 


100 


500 


Elliptical (bipolar core?)'' 


LMC-SMP 83 


95 


170 


1600 


Spherical*, Bipolar Ring* 


LMC-SMP 29 




200 






SMC-SMP 1 


37 


41 


500 


Spherical 6 '', Round" 


SMC-SMP 3 


95 


92 


> 2000 


Bipolar Ring 6 '', Bipolar w/bipolar core" 


SMC-SMP 22 




200 




Bipolar Ring 6 , Bipolar w/bipolar core" 



REFERENCES — (aa): This p a per, (bb): Dopita & Meatheringhanj I1991al lbh: 
Dopitaet al.1 f l997>: Vassiliadis et al. (1996, 1998a) (a): Stanghell jnTet1u'Ml999|) (b): 
Vassihadis et jriT 7l998b) (c): D opita etalH 19961) (dl: IShaw et al.l 4200lf) 



TABLE 8 

Levels and superlevels for model ions 



Elemen t 


I 


II 


in 


IV 


V 


VI 


VII 


VIII 


IX 


X 


XI 


H 


20/30 


1/1 




















He 


40/45 


22/30 


1/1 


















C 






30/54 


13/18 


1/1 














N 








29/53 


13/21 


1/1 












O 








29/48 


41/78 


13/19 


1/1 










Si 








22/33 


1/1 














p 








36/178 


16/62 


1/1 












S 








51/142 


31/98 


28/58 


1/1 










Fe 








51/294 


47/191 


44/433 


41/254 


53/324 


52/490 


43/210 


1/1 



